Gate tunable infrared phonon anomalies in bilayer graphene 



0^ 
O 
O 
(N 

00 



I 

o 
o 



> 
m 
o 

(N 

o 

On 
O 



A. B. Kuzmenko,^ L. Benfatto,^'^'^ E. Cappelluti,'^' I. Crassee/ 
D. van der Marel,^ P. Blake,^ K. S. Novoselov,^ and A. K. Gcim^ 

^ Departement de Physique de la Matiere Condensee, 
Universtte de Geneve, CH-1211 Geneve 4, Switzerland 
^Centra Studi e Ricerche "Enrico Fermi", via Panispema 89/A, 1-00184, Rome, Raly 
■''SMC Research Center, CNR-INFM, and ISC-CNR, via dei Taurmi 19, 00185 Roma, Italy 
^Dipart. di Fisica, Universita "La Sapienza", P.le A. Mora 2, 00185 Rome, Raly 
^Manchester Centre for Mesoscience and Nanotechnology, 
University of Manchester, Manchester MIS 9PL, UK 

We observe a giant increase of the infrared intensity and a softening of the in-plane antisymmetric 
phonon mode E„ (~ 0.2 eV) in bilayer graphene as a function of the gate-induced doping. The 
phonon peak has a pronounced Fano-hke asymmetry. We suggest that the intensity growth and 
the softening originate from the coupHng of the phonon mode to the narrow electronic transition 
between parallel bands of the same character, while the asymmetry is due to the interaction with 
the continuum of transitions between the lowest hole and electron bands. The growth of the peak 
can be interpreted as a "charged- phonon" effect observed previously in organic chain conductors 
and doped fuUerenes, which can be tuned in graphene with the gate voltage. 



There is a lot of interest in electronic properties of 
graphene, especially, due to the chiral nature of its quasi- 
particles pj]. Graphene also allows tuning the Fermi 
level continuously through zero energy by the ambipo- 
lar electric field effect, and this tunability has been in- 
tensively exploited in many experiments. In particular, 
it was found that electric field strongly affects electron- 
phonon interactions [1, 0, 0, S H H [o| , which is hard 
or impossible to achieve in bulk materials. Particularly 
intriguing is the case of bilayer graphene, where the elec- 
tric field breaks down the interlayer symmetry f?| and 
changes qualitatively the electron-phonon coupling [9]. 
Up to now the vast majority of the optical phonon stud- 
ies in graphene were performed by Raman spectroscopy, 
whereas infrared techniques that can provide important 
complementary information, due to different selection 
rules, remain largely unexplored. Here we present an 
infrared study of the in-plane optical phonon mode in 
gated bilayer graphene. We observe several new phonon 
anomalies, of which the most spectacular one is a giant 
enhancement of the phonon resonance as a function of 
the gate voltage. In addition, the phonon peak is found 
to have a pronounced Fano shape 10], which indicates 
a coupling of this mode to a continuum of electron-hole 
excitations. We believe that this gate-tunable coupling 
can be used in various electro-optical applications. 

Reflectivity spectra of an exfoliated bottom-gated bi- 
layer graphene flake on top of a SiO2(300 nm)/Si sub- 
strate (Fig. la) were taken at 10 K at near-normal inci- 
dence using an infrared microscope attached to a Fourier 
transform spectrometer. Lithographically made leads al- 
lowed applying the gate voltage Vg and measuring resis- 
tivity during the optical experiment. The charge neutral 
state corresponding to the resistivity maximum was ob- 
served at Vgo K, -30 V. Fig. 2a shows the reflectivity 
normalized to the bare substrate at various gate voltages 
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FIG. 1: (Color online) (a) A sketch of the infrared reflectivity 
experiment and the crystal structure of Bernal-stacked bilayer 
graphene. Atomic displacements corresponding to the anti- 
symmetric mode Eu are indicated by arrows; for the symmet- 
ric mode Eg the displacements in one of the layers should be 
inverted, (b) A simplified band structure of bilayer graphene. 
Arrows indicate electronic transitions discussed in the text. 



from -100 V to 100 V. The doping level can be deduced 
using the relation n[cm.-'^] = 7.2-10i° x(Vg - Vgo) [V]. 
One can clearly see a phonon structure at about 1600 
cm^^ (0.2 eV) corresponding to the infrared active anti- 
symmetric Eu mode (Fig. la). In order to easier compare 
the phonon feature at different values of Vg, in Fig. 2b 
we subtract from the reflectivity curves a smooth base- 
line obtained by their low-order polynomial fltting out- 
side the region 1500 - 1700 cm^^. While the structure 
is weak in the vicinity of Vg = Vgo, it grows quickly at 
elevated (both positive and negative) gate voltages. 

It is not straightforward to interpret the phonon line- 
shapes in the raw reflectivity spectra, which are sensi- 
tive to the choice of the substrate. In Fig. 3 we show the 
real part of the optical conductivity of graphene, A(Ti(w), 
extracted using a Kramers-Kronig consistent variational 
spectra analysis llj that takes into account the known 
optical constants and thicknesses of the Si02 and Si lay- 
ers. The smooth background due to direct optical ab- 
sorption by electronic transitions, which is discussed in 
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FIG. 2: (Color online) (a) infrared reflectivity spectra of bi- 
layer graphene on SiO2(300 nm)/Si around the Eu peak nor- 
malized by the reflectivity of the bare substrate at different 
gate voltages, (b) the same spectra with a baseline sub- 
tracted. The curves (except at -fOO V) are vertically shifted. 



detail elsewhere [l^ll3|, is subtracted. 

A few interesting observations can be made. First of 
all, the peak intensity grows rapidly as the gate voltage 
is tuned away from Vgo. Second, the phonon peak is 
clearly asymmetric, with a characteristic Fano lineshape 
[lot that is most pronounced at not very high values of 
\Vg-Vgo\ (for example, at Vg = -70, -80, 0, +10 V). Third, 
close to the charge neutral point (for example at -40 and 
-50 V), the peak is negative, i.e. the conductivity has a 
weak but a distinct dip. Fourth, the position of the peak 
maximum is gate-voltage dependent. For example, one 
can see that for Vg > 40 V it shifts to lower frequencies 
as the gate voltage is increased. 

The Fano lineshape provides evidence of a coupling 
between this sharp mode and a continuum of electronic 
excitations in the same energy range. Following the gen- 
eral Fano theory the optical conductivity Aai(w) of 
such coupled system minus the background due to the 
bare continuum can described by the formula: 



Acri(u;) 



2qz - 1 



47rr q^{l + z^) 



(1) 



with z = 2{uj—ujq) /F. Here is the phonon frequency, F 
is the linewidth and LOp is the plasma frequency, which is 
proportional to the infrared effective charge as specified 
below. The asymmetry is described by a dimensionless 
parameter q so that the symmetric Lorenzian lineshape is 
recovered in the limit \q\ — > oo. The doping dependence 
of the Fano parameters extracted by the fitting of the 
conductivity spectra of Fig. 3 is shown in Fig.4a-d. 

From Fig. 4a one can see that ojq decreases almost lin- 
early with the slope divo/dVg « -0.058 cm~^/V. Note 
that due to the Fano shape, ujq is larger than the posi- 
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FIG. 3; (Color online) (a) The real part of the optical conduc- 
tivity Acri(aj) of bilayer graphene at different gate voltages. 
The electronic baseline is subtracted. The curves (except at - 
100 V) are vertically shifted, (b) and (c): the expanded views 
of the Eu peak at electron and hole doping. 



tion of the conductivity peak maximum, the difference 
being dependent on the asymmetry parameter. This is 
why the softening is seemingly less pronounced in FiglSlD 
and [St. The linewidth is enhanced at low doping levels 
and decreases with the increasing the gate voltage. The 
fitted value of F shows a saturation at 10 cm~^, which 
is the resolution limit of this measurement and there- 
fore an overestimation of the true linewidth. The plasma 
frequency ujp is at minimum close to Vgo (Fig. 4c), where 
ujp ^ 200 cm"-'^ and strongly increases with the gate volt- 
age. In Fig.4e we plot the bare phonon spectral weight 
W = ujp/8 as well as the peak area W = J Aai{ui)du} 
(integrated between 1570 and 1610 cm^^). Note that 
in the presence of coupling W « — l/g^)/8 ^ W. 
If the asymmetry is very strong {\q\ < 1) then the area 
becomes negative. This is because W is affected not 
only by the direct optical absorption due to the phonon 
mode but also by the coupling-induced transformations 
of the optical response of the continuum. In particular 
the negative value of W is explained by the removal of 
spectral weight from the continuum. The asymmetry pa- 
rameter |g| also shows a strong doping dependence, with 
a decrease as the doping is increased. This decrease is 
not monotonic: the strongest changes are observed at 



Vn-Vr 



50 V, while at higher doping the asymmetry 



remains roughly constant. 
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FIG. 4: (Color online) (a)-(d) Fano parameters loo, F, uip (^off) 
and q of the phonon peak as a function of doping. At very 
low doping, the error bars for ujq exceed the panel scale. The 
effective charge corresponding to the calculated static dipole 
is shown by the dashed line in panel (c). (e) The bare spectral 
phonon weight W and the peak area W' . (f) The background 
electronic conductivity at the phonon energy. 



Let us discuss the observed phonon anomalies and their 
relation to the band structure and electronic excitations. 
The frequency shift Awq and the linewidth F are related 
to the phonon self energy n(a;): Aujq « h~^Re n(wo), 
r « —h~^lm n(wo). The self energy for the zone-center 
optical phonons is given at zero temperature by a sum of 
the contributions of all vertical (momentum conserving) 
transitions across the Fermi level weighted with their re- 
spective electron-phonon matrix elements. Each contri- 
bution has a Lorentzian spectral shape (Ae — hui + iS)~^, 
where Ae is the transition energy and 5 is a broaden- 
ing parameter. Therefore F is determined only by the 
transitions that overlap with the phonon energy within 
6, while Awo is affected by all transitions with non-zero 
electron-phonon matrix elements. 

A simplified band structure of bilayer graphene, where 
we neglect the electron-hole asymmetry and the bandgap, 
is shown in Fig. lb. Each of the hole and electron bands 
is split into two parallel subbands (numbered 1,2 and 3,4 
respectively) separated by the interplane hopping energy 
71 ~ 0.4 eV. We shall refer to the case of electron doping 
(positive gate voltages) although the results are equally 
valid for hole doping. Based on this band structure, Ando 
calculated the phonon self energy and predicted the soft- 
ening of the phonon with doping f7^, in an excel- 



lent agreement with the present infrared measurement. 
Interestingly, the softening was also found in a Raman 
study Q, where the mode was observable due to a 
possible breaking of the inversion symmetry of bilayer 
graphene in the used device. By comparing the observed 
slope dujo I dVg with the theoretical curve we estimate the 
dimensionless electron-phonon coupling parameter A de- 
fined in RefQ to be about 6 x 10 ^, in agreement with 
the Raman study [1] . According to Ando [TI] , the soften- 
ing of the Eu mode is due to the contribution to the self 
energy coming from the transitions between the bands 3 
and 4 (thick red arrows in Fig. lb). Since the two bands 
are parallel, these transitions form a narrow peak cen- 
tered at 7i. Its intensity grows rapidly with doping in 
the course of the progressive filling of the band 3. 

It is important that none of the interband transitions, 
except the ones between the bands 2 and 3, overlap with 
the phonon energy, as 71 is about twice as large as cjq. 
On the other hand, the electron-phonon matrix elements 
between the £"„ mode and the 2 — > 3 excitations van- 
ish according to the symmetry analysis 0]. Therefore, 
the same calculation [7| predicts a very small and weakly 
doping dependent linewidth of the E^ phonon mode. In- 
stead, we observe that F is large (~ 20-30 cm"-*^) at low 
doping and strongly decreases with n (Fig. 4b). 

We argue that the large linewidth and the Fano shape 
of the phonon peak are likely due to the interaction with 
the 2 — > 3 transitions (shown with thin green arrows in 
Fig. lb), despite the fact that the self-energy calculation 
Ref.01 would suggest a negligible coupling in this case. 

The 2^3 transitions form a continuum of particle- 
hole excitations. Such a continuum is also present in 
monolayer graphene and in graphite, where it results 
in the universal optical conductivity of e^/Ah [l 
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IT^ [l7j . The continuum starts at u; = 2Ep due to the 
Pauli principle [l3|. If Ep < wo/2 then the phonon 
frequency lies inside the 2 — > 3 continuum and in this 
case one expects the linewidth and the asymmetry to be 
stronger. Indeed, we see a clear correlation between the 
sudden decrease of the asymmetry parameter q (Fig.4d), 
the decrease of the broadening (Fig. 4b) and the reduc- 
tion of the background electronic optical conductivity at 
the phonon frequency abg{ijJo) (Fig-4e) [l^], which corre- 
sponds to Ep K, ojq/2. It is worth mentioning that the 
anomalous softening of the Raman-active phonon mode 
in bilayer graphene is observed at about the same gate 
voltage 0]. 

Among the possible reasons of the " forbidden" interac- 
tion between the mode and the 2^3 electronic exci- 
tations we can mention (i) the electron-hole asymmetry, 
which is clearly observed in the electronic infrared spec- 
tra [ll[ll[ll[ll] and (ii) an electric field generated by 
the gate electrode and charged impurities which breaks 
the inversion symmetry Note that a Fano phonon 
lineshape in the infrared spectra of bilayer graphene was 
also observed in Ref . |20l . 
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Finally, we discuss the most spectacular effect of this 
study: a large increase of the phonon intensity as a func- 
tion of doping. From the plasma frequency we can de- 
duce the dimensionless effective infrared charge of the 
mode Zcs = Wp(mcT^/(47r7Ve^))^/^, where = 4 is the 
number of atoms in the unit cell, V is the unit cell vol- 
ume, nic is the mass of the carbon atom and e is the 
elementary charge. The doping dependence of the effec- 
tive charge is shown in Fig. 4c. One can see that at high 
gate voltages Z^e exceeds 1, which is the value observed 
in entirely ionic compounds like NaCl and LiF. We first 
compare this value to the effective charge of the static 
dipole created because the doped charge is distributed 
not equally between the atomic positions A1(A2) and 
B1(B2) (Fig. la). The static charge calculated using the 
tight-binding model (the dashed line in Fig. 4c) appears 
to be more than 300 times smaller than the observed 
infrared (dynamical) effective charge. Such an utter dis- 
crepancy calls for a radically different physical approach. 

Notably, similar effects of a large enhancement of vi- 
brational resonances were observed in organic linear- 
chain conductors such as K-TCNQ [2l| and in doped Ceo 
[2^ . Effective infrared charges significantly larger than 
ionic charges were also found in some transition metal 
compounds, for example in FeSi (23l |. These findings were 
explained using the so-called " charged-phonon" theory of 
M.J.Rice et al. 2j|. Within this model, the molecular 



vibrations that are usually not active or weakly active in 
the infrared range become enhanced due to a coupling to 
electronic resonances at higher energy (such as the tran- 
sition tiu tig in doped Cgo)- Their effective charge is 
effectively "borrowed" from the electronic transitions, in 
accordance with the f-sum rule. 

A clear correlation between the intensity of the phonon 
resonance and the strength of the electronic transition 
was experimentally observed 22|. The charged-phonon 



model thus seems to be relevant to the intensity enhance- 
ment of the Eu mode in bilayer graphene, in view of the 
observation that the strong doping dependence of the in- 
tensity of the 3^4 transition correlates with the growth 
of the phonon peak. Interestingly, a simultaneous in- 
crease of the phonon intensity and electronic absorption 
with doping was also found in organic field effect tran- 



sistors [25|- We note that even at zero doping, where 
the intensity of the 3^4 transition is zero, the effective 
charge is about 0.25. This suggests that other transitions 
(1 — > 3, 2 — >• 4 and 1 — > 4), which remain populated at 
all doping levels, also contribute to Zee- We notice that 
the zero-doping value of Z^s is comparable with the ones 
reported in graphite [2^ . 

As a concluding remark, the possibility to control the 
phonon intensity by the gate voltage might be of practi- 
cal importance, for example, in optoelectronics. We also 
suggest that the peak can be used as a doping indicator. 

In summary, we observed several anomalies of the in- 
frared active mode in bilayer graphene: (i) a drastic 



intensity enhancement and a softening as a function of 
doping, (ii) a strong Fano-like asymmetry of the line- 
shape, with even a negative peak area at low doping. We 
suggest that the coupling to the narrow transition 3^4 
is responsible for the mode softening ^ and the inten- 
sity enhancement, via the charged-phonon mechanism 
241 . The interaction with the continuum 2 — » 3 likely 
gives rise to the asymmetry and the broadening of the 
phonon peak. Further studies are required to understand 
quantitatively the complicated picture of the electron- 
phonon coupling in graphene. This work was supported 
by the Swiss National Science Foundation (SNSF) by 
the grant 200021-120347, through the National Center 
of Competence in Research "Materials with Novel Elec- 
tronic Properties-MaNEP" and by Italian MIUR project 
PRIN 2007FW3MJX. 
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